NASA TECHNICAL NOTE 



1 NASA TN D-5460 


LOAN COPY: RETURN TO 
AFWL (WlA-2) 
KIRTLANO AFB, N MEX 



INJECTION OF AN INVISCID 
SEPARATED JET AT AN OBLIQUE 
ANGLE TO A MOVING STREAM 

by Marvin E. Goldstein and Willis Braun 

Lewis Research Center 
Cleveland, Ohio 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • OCTOBER 1969 


tech library kafb, nm 



TECH LIBRARY KAFB, NM 


1 . Report No. 

NASA TN D-5460 

4. Title and Subtitle 


2. Government Accession No. 


INJECTION OF AN INVLSCID SEPARATED JET AT AN 
OBLIQUE ANGLE TO A MOVING STREAM 



3. Recipient's Catalog No. 


5. Report Date 

October 1969 

6. Performing Organization Code 


7. Author(s) 

Marvin E. Goldstein and Willis Braun 

9. Performing Organization Name and Address 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 44135 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 

15. Supplementary Notes 


16. Abstract 


8. Performing Organization Report No. 

E-5039 

10. Work Unit No. 

129-01-07-07-22 

1 1 . Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Note 


14. Sponsoring Agency Code 


An analytical solution has been obtained to the problem of a two-dimensional incompress- 
ible jet injected into a moving stream from an orifice set at an oblique angle to the 
stream. It is assumed that the jet separates from the downstream edge of the orifice. 
The solution is valid when the normalized difference between the total pressure in the 
jet and the total pressure in the main stream is not too large. Typical flow patterns are 
shown to illustrate the effects of varying both the jet offset ratio and the total pressure 
within the jet. 


17. Key Words (Suggested by Author(s)) 


Jet injection 
Inviscid flow 
Complex variable 
Conformal mapping 


Unclassified 

- unlimited 

19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

Unclassified 

Unclassified 

71 


18. Distribution Statement 


22. Price* 

$3.00 


*For sale by the Clearinghouse for Federal Scientific and Technical Information 

Springfield, Virginia 22151 



INJECTION OF AN INVISCID SEPARATED JET AT AN OBLIQUE 
ANGLE TO A MOVING STREAM 
by Marvin E. Goldstein and Willis Braun 
Lewis Research Center 

SUMMARY 

An analytical solution has been obtained to the problem of a two-dimensional incom- 
pressible jet injected into a moving stream from an orifice set at an oblique angle to the 
stream. It is assumed that the jet separates from the downstream edge of the orifice. 
The solution is valid when the normalized difference between the total pressure in the 
jet and the total pressure in the main stream is not too large. Typical flow patterns are 
shown to illustrate the effects of varying both the jet offset ratio and the total pressure 
within the jet. The analysis shows that for orifices tilted into the main stream small in- 
creases in the total pressure in the jet result in large increases in the jet penetration 
and thickness. 


INTRODUCTION 

The flow field which results from the oblique injection of a jet into a moving stream 
is of considerable interest in a number of fluid mechanical devices. Among these are 
ground- effect machines, jet flaps, wing fans on VTOL aircraft and fuel injection systems 
in combustion chambers. 

The fluid mechanics of jet penetration into moving streams is by no means fully 
understood. Some insight into this phenomenon can be gained by considering the injec- 
tion of two-dimensional inviscid incompressible jets into moving streams since such 
flows are simple enough to be amendable to mathematical analysis. Of course, viscous 
effects can be quite significant in real fluid flows. However, in order to take viscous 
effects into account the usual procedure is to first perform an inviscid analysis and then, 
provided viscous effects do not modify the flow significantly in the wake region, to mod- 
ify the flow by superimposing viscous boundary layers. In any event, it is hoped that the 
inviscid analysis will reveal some of the important features of the flow field and thereby 



lead to an increased understanding of the phenomena involved. 

Relatively few analyses of inviscid flows of this type have been performed. This is 
due at least partially to the fact that analysis of such flows involves the solution of a 
nonlinear problem in which the shapes of the boundaries are unknown. However, in the 
very special case where the total pressure in the jet is equal to the total pressure in the 
main stream, the classical theory of inviscid jets can be used to obtain solutions. Flows 
in which the total pressure in the jet equals the total pressure in the main stream are 
discussed by Ehrich in reference 1. Ehrich considers a large number of possible flow 
configurations for jets issuing from both slots and orifices. It has been found, however, 
that there is a serious error in his results for jets issuing from orifices. (Since these 
results emerge as a special case of the analysis performed herein the corrected solu- 
tions will be given.) The error was caused by taking the wrong sign for the square root 
in the expression for the contraction ratio. 

There is another limiting case, called the ” strong jet approximation, " in which the 
analysis can be considerably simplified. This is the case where the total pressure in 
the jet is very much larger than that in the main stream. An analysis of this type of jet 
was first carried out by G. I. Taylor (ref. 2). Taylor obtained an analytical solution by 
introducing an additional approximation. This latter approximation was removed by 
Ackerberg and Pal (ref. 3), who developed a variational principle for the problem and 
thereby obtained a numerical solution. Although making the strong jet approximation 
considerably simplifies the analysis, the resulting problem is still nonlinear and in- 
volves an unknown boundary so that an exact analytical solution does not seem feasible. 

In any event, as was pointed out by Taylor (ref. 2), the viscous spreading of real jets 
which have a total pressure much larger than that in the main stream is so large that the 
inviscid solutions show no relation to actual experiments. 

The general problem of an inviscid jet issuing into a flowing stream from a two- 
dimensional vaned slot was reformulated by Ting, Libby, and Ruger (ref. 4) in terms of 
two simultaneous nonlinear singular integral equations. Because of the extreme diffi- 
culty involved in solving such equations, the authors considered two limiting cases of the 
equations. The first corresponded to the strong jet approximation and the second to the 
case where the total pressure in the main stream differs by only a small amount from 
the total pressure in the jet. Because of the nature of the boundary conditions associ- 
ated with the vaned slot, the authors could not linearize the problem even in the case of 
small total pressure difference. Hence, their formulation of this problem is still in 
terms of a very difficult nonlinear singular integral equation. The authors give a nu- 
merical procedure for obtaining a solution to this equation; but, as pointed out by Ting 
and Ruger (ref. 5), no attempts at carrying out the solution have been successful. It 
was also shown in reference 5 that no ordinary perturbation procedure could be used to 
obtain the solution to the problem of a jet issuing from a vaned slot. However, it will 
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be shown subsequently that a systematic perturbation procedure can be used to linearize 
the problem of a jet issuing from an orifice into a moving stream in the case where the 
normalized total pressure difference is small. 

In this report an explicit analytic (closed form) solution is obtained for the flow field 
resulting from a jet issuing from an orifice into a moving stream. The orifice is set at 
an oblique angle to the flow and it is assumed that the jet separates from the downstream 
edge of the orifice to form a stagnent wake. The flow configuration is shown in figure 1. 
The flow will be assumed to be two-dimensional, inviscid, and incompressible. In addi- 
tion, it will be required that in a certain sense (to be specified more precisely below) 
the normalized difference between the total pressure in the jet and the total pressure in 
the main stream be small. The upstream boundary of the jet is the stream line emanat- 
ing from the upstream edge of the orifice. The velocity, in general, will not be continu- 
ous across this stream line. 

The problem is solved by expanding the solutions in a small parameter related to 
the difference in total pressure between the jet and the main stream. The zeroth-order 
solution corresponds to equal total pressures and can be written down immediately by a 
simple application of classical techniques. The solution to the zeroth-order problem in 
fact corresponds to one of the solutions obtained in reference 1. However, it was found 
that by using a somewhat different procedure than that used in reference 1 a simpler and 
more convenient form of this solution could be obtained. 

Since the boundary shapes for the first-order (different total pressure) problem are 
unknown, a technique similar to that employed in thin airfoil theory is used to transform 
the first-order boundary conditions to the zeroth-order boundary. This transformed 
problem still involves a combination of boundary and jump conditions which cannot be 
handled by ordinary techniques. Therefore, a new procedure was developed to trans- 
form this problem into a standard problem for a sectionally analytic function. (A sec- 
tionally analytic function is one which is holomorphic in each of two adjoining regions 
and has a specified jump in value across the boundary of these regions. ) The procedure 
consists of introducing a new dependent variable in such a way that the new variable has 
to satisfy only jump and symmetry conditions instead of the combination of jump condi- 
tions and boundary conditions that the original variable satisfies, hi order to introduce 
this new variable, several mappings between certain complex planes are introduced. 

The solution is then obtained by using the theory of sectionally analytic functions. 

It should be emphasized that the techniques developed herein are quite general and 
can be applied to a wide variety of jet injection problems. Since it is impossible to tell 
from the inviscid analysis whether separation will occur at the downstream edge of the 
orifice, the case with no separation will be considered in a future report. It is shown in 
appendix A that, by a simple rescaling, the results obtained herein can be applied to the 
case where the density of the fluid in the jet is different from that in the main stream. 
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It is perhaps worth pointing out why an ordinary perturbation procedure works for 
jet injection problems from orifices but will not work for jet injection problems from 
slots and vaned slots. The reason is a consequence of the fact that in the mathematical 
solution the flow fields associated with slots and vaned slots have stagnation points at the 
up-stream edge of the slot, whereas the flow fields associated with orifices do not. The 
nature of this stagnation point changes discontinuously as the total pressure in the jet is 
changed from the free-stream value. This discontinuous change results in a singularity 
at the stagnation point in the first-order solution when an ordinary expansion procedure 
is used. This difficulty can be overcome by using the method of matched asymptotic ex- 
pansions. The inner expansion can, in fact, be obtained by rescaling the results given 
in reference 5, and the outer expansion can be obtained by the procedure given in this 
report. 


ANALYSIS 

Formulation and Boundary Conditions 

It will be assumed that the flow is inviscid, incompressible, and irrotational. The 
jet configuration is illustrated in figure 1. The analysis is limited to the case in which 
the difference between the total pressure in the jet Pj and the total pressure in the main 
stream is not too large; or more specifically, to the case in which 



Figure 1. - Schematic configuration of jet penetrating into flowing stream. 
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is the density of the fluid, and is the velocity of the main stream at infinity. 

Let l be a convenient reference length which will be specified in the course of the 
alysis. The X and Y components of the velocity, U and V respectively, will be 
made dimensionless by V^, and the stream function ^ and the velocity potential $ 
will be made dimensionless by V^l. Thus, the dimensionless quantities u, v, and 
<p are defined by u s U/V^, v 5 V/V M , & = and cp = #/V M Z. The dimension- 

less complex conjugate velocity £ and the dimensionless complex potential W are de- 
fined, as usual, by 


£ = u - iv 


and 


W - (p + 


With all lengths made dimensionless by l (i.e., x = X/l, y = Y/l, a = A/l, b = B/l, 
and h = B/l), the flow configuration is shown in the physical plane (with the complex 
variable z defined by z = x + iy) in figure 2. 

The stream of fluid issuing from the orifice formed by the two parallel walls KD 
and EK meets the main stream at the point D and forms a common stream line, which 
is denoted by S in figure 2. At the point E the jet separates from the wall EK and a 
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stream line (denoted by F in fig. 2) marks the boundary between the jet and the down- 
stream wake region, which is assumed to be composed of stagnent fluid. The shapes of 
the curves S and F are a priori unknown and must be determined from the analysis. 

s s s 

Points on the curves S and F will be denoted, respectively, by z = x + iy and 

F F . F 
z + x + iy . 

In order to satisfy the requirement that there be no discontinuities in static pressure 
anywhere within the flow field, it is necessary (as will be shown subsequently) to allow 
the velocity to be discontinuous across S. For this reason the stream line S will be 
called the slip line. The region within the jet and orifice is denoted by D + , and the re- 
maining region of the flow (i. e. , the main stream) by D~. Since the velocity (and as a 
consequence, the velocity potential) is discontinuous across S, it is convenient to use 
the superscripts + and - to distinguish these regions. Thus, 


S(z) = 


C + (z); z e D 4 


C(z); z e D" 


and 


W(z) = 1 


W + (z); z G D 4 


W"(z); z G D" 


Then and W + are holomorphic in the interior of D + , and and W~ are holo- 
morphic in the interior of D”. 

The following argument will show that the velocity must be discontinuous across S. 
Bernoulli's equation can be written for the flow inside the jet as 


P(x, y) 

Vi 


+ k + (z)i 2 


p i 


L pvl 


D 


( 2 ) 


and for the flow external to the jet, Bernoulli's equation is 


m + |r(z)i 2 = 


D 


ipV* 

2 


-pvl 

2 


( 3 ) 
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Now, in view of the fact that S is a common stream line to the external and internal 

g 

flow, the points z in both these equations can be chosen to be the same point z on the 
curve S. Since the static pressure must not be discontinuous across S, it follows from 
equations (2) and (3) that 


|;V>| 2 - | c V s ) | 2 



(4) 


g 

at every point z of S. 

On the other hand, since there is no flow in the wake, it is clear that the static pres- 
sure along the free- stream line F is constant and equal to the pressure in the wake. 
Because the velocity in the jet must become uniform across the jet far downstream (i. e. , 
at the point C in fig. 2), it follows that there cannot be any pressure difference across 
the jet far downstream. It can be concluded from this that the pressure in the wake must 
equal p , the static pressure of the main stream at infinity. Hence, equation (2) shows 
that at any point z of F 


l? + (z F )| 2 





= 1 + 6 


- pV 2 
2 P °° 


(5) 


Since S is a common stream line to the internal and external flows it is clear that 
y*tiW*(z S ) and y#tW-{z S ) are both constants. Moreover, the arbitrariness in the defini- 
tion of W can be partly removed by choosing these constants to be zero (ref. 6). Hence, 

^y»W + (z) = (z) = 0 for z e S (6) 

The remaining arbitrariness in W can be removed by choosing 

W+(0) = W"(0) = 0 (7) 

F 

The fact that F is a stream line implies that there exists a real constant if/ such that 

j?,»W + (z) = if/ F for z e F (8) 

The conditions imposed on the velocity at infinity are (in view of the manner of nondimen- 
sionalization; also see fig. 2) 
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£ + (z) - 0 for z — K 


( 9 ) 


and 


£ - (z) — 1 for z — G 


( 10 ) 


The remaining boundary condition is that the normal component of the velocity vanish on 
the solid boundaries. These conditions are sufficient to completely determine the solu- 
tion. They are summarized below for convenient reference. 


|? + (z)| 2 - |r(z)| 2 = e "I 



r z 

9mW ¥ {z) = 9mW~(z) 

' = °J 

|? + (z)| 2 = 1 + 

>z e F 

9»tW + (z) = i// F 

J 


9M''C(z) = 0; 

z e KD 

9mC(z) = 0; 

z e EK 

jWr(z)=0; : 

zGGD 


( 11 ) 


Asymptotic Expansions 

For small values of e, the functions £* and W* can be expanded in an asymptotic 
power series in e. In view of the fact that the shape of the slip line and of the free- 
stream line depend on e, these expansions imply that the coordinates of S and F, z® 
and z , respectively, and the asymptotic jet width h must also be expanded in powers 
of e. Hence, 
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o 

II 

+ +. 

II 

o 

+ eW* + * 

s 

S 

~ z o 

+ ez 1 + , 

F 

F , 

= z o 

+ ezj + . 

O 

■d 

II 

+ ehj + . 

fe o 

II 

F 

+ + . 

expansions for 


( 12 ) 


complex variable z is being expanded. These expansions are only for the shape of the 
boundary of the jet. In fact, only one of the expansions (12) is independent, and the 
remaining expansions are determined in terms of it. Thus, for example, since 


t* * m 


dz 


it follows that 


£q + ^ + • • 


dW, 


dz 


+ e 


dW, 

dz 


Hence, equating like powers of e shows 




dz 


,± ^ 


• • • 


(13) 


I 



Therefore, once the coefficients of the various powers of e are known in the expansion 
for £*, the coefficients of the various powers of e in the expansion for W ± are deter- 
mined by equations (13). 

The reason for omitting the superscript + or - in the zeroth-order terms of the 
first two expansions is that (as will be shown subsequently) the zeroth-order solutions 
are not discontinuous across the curve S and so there is a single function £q which is 
homomorphic in the entire flow field (of course, the same is true for Wq). 

The reference length Z will now be chosen in such a way that 


Thus, Z is the zeroth-order asymptotic thickness of the jet. This is denoted symbolic- 
ally by putting 


Z 



(14) 


The last expansion (eq. (12)) is then 

h = 1 + eh^ + . . . (15) 

Zeroth-Order Solution 

When the expansions (12) are substituted into the boun da i y conditions (7) 
and (9) to (11)) and only the zeroth-order terms are retained, the following boundary 
conditions for the zeroth-order solution are obtained: First, the first boundary condi- 
tion (11) shows, as has already been anticipated, that the zeroth-order solution 
must be continuous across the slip line and, hence, that it is characterized by functions 
which are holomorphic everywhere within the flow field. The remaining conditions show 
that 


W Q (0) = 0 





(16) 
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| £ o( z o)| 

-1 

^ w o( z o 

II 

o%, 

V 


'z e KD 

9*nt 0 ( z ) = o< 

z e EK 


^z ^ GD 


£ 0 (z) - 0 for z - K 

£ 0 (z) - 1 for z — 


( 17 ) 


The conditions (16) merely serve to show that, because of the manner in which the 
arbitrary constants have been adjusted in the complex potential, the stream line ema- 
nating from the point D is to be taken as the zero stream line. 

Now the change in the stream function across the jet must be equal to the volume 
flow rate through the jet. Hence, if Qq denotes the zeroth-order volume flow through 
the jet, in view of the second condition (16) and the second condition (17), it 
follows that 


% = 0 - *0 ■ 


The first boundary condition (eq. 17)) shows that far downstream in the jet (i. e. , at the 
point C) the zeroth-order velocity goes to 1. In view of the normalization (15) the as- 
ymptotic thickness of this portion of the jet must also be 1. It follows from these re- 
marks that Q q = 1. Hence, 

*0=-l ( 18 ) 

Now, the boundary value problem posed by the boundary conditions (17) is a simple 
free-stream-line problem which can be readily solved by the Helmholtz- Kir chhoff tech- 
nique. In fact, the solution to this problem has already been carried out by Ehrich 
(ref. 1). His solution, however, is somewhat inconvenient for our purposes. The pro- 
cedure used herein for obtaining the solution is (ref. 6) to draw the region of flow in the 
hodograph plane and in the complex potential plane, and then to find the appropriate map- 
ping of these two planes into some convenient intermediate plane (say the T-plane). The 
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Figure 3. - Zeroth-order complex potential plane (Wg-plane). 



shapes of these regions can readily be deduced from the boundary conditions (17), and 
they are shown in figures 3 and 4 (we have put Wq = <Pq + i’/'Q in fig. 3). The corre- 
sponding points in the various planes are designated by the same letters. The zeroth- 
order "slip line" is shown dashed in these figures since it does not correspond to a line 
of discontinuity and can therefore be ignored as far as obtaining the zeroth-order solu- 
tion is concerned. The intermediate T-plane is chosen in such a way that the region of 
flow maps into the upper half-plane in the manner indicated in figure 5. We shall denote 
the real and imaginary parts of the variable T by ^ and rj, respectively. The region 
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Figure 5. - Intermediate plane (T-plane). 


of the T-plane into which the zeroth-order flow field interior to the jet and orifice maps 
is denoted by and the region of the T-plane into which the zeroth-order main 
stream maps is denoted by q. The dividing line between these two regions (which is 
being called for convenience the zeroth-order slip line even though no slip occurs in the 
zeroth-order solution) is denoted by y^. 

A simple application of the Schwartz -Chris toff el transformation (ref. 7) shows that 
the mappings which properly transform the Wq- plane and the fQ-plane into the upper 
half T-plane in the manner indicated in the figures are respectively defined by 


and 


dWQ 

dT 


1 T + 1 

7 T T 


?7 > 0 


(19) 


d In 
dT 


iVX 

TVT - A 


7] — 0 


Or, performing the indicated integrations, 

W 0 = - (T + 1 + In T) - i 77 > 0 

U IT 


( 20 ) 


( 21 ) 
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rj > 0 


( 22 ) 


_ T - 2 A - 2iVZ Vt - A 


1 _ T - 2 A + 2iVA Vt - A 

^o" T 


In order to find the equation for SSq, the zeroth- order slip line in the T- plane, notice 
figures 3 and 5 show that ,/*#Wq(T) = 0, whenever T e y^. Now equation (21) shows, 
since T = £, + it], 


W q (T) = i 


| + 1 + it] - in + In ^ + 77^ 


1/2 


+ i tan" ^ 

L 


0 < tan" < 7r 


Hence, setting the imaginary part of this expression equal to zero we find that 


1 

77 


(” 


77 + tan" * 


0 


for T e y Q 


This equation has two solutions. The first 

77 = 0 and | < - 1 


corresponds to the negative real axis in the Wg-plane and is of no interest to us. The 
other solution is 


tan ( 77 -77) =2 0 < 77 < 77 

I 

This solution corresponds to the positive real axis in the Wg-plane and therefore deter- 
mines the equation of y q. Thus, 

£ = U = 0 < 77 < 77 for TGy fl 

tan(77 - 77) tan 77 

Hence, the parametric equation for y^ is 

T = - + xt] = - -J— e" i7? 0 < 77 < 77 (24) 

tan 77 sin 77 
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It will be necessary in what follows to have an expression for VT - A for T e in 
terms of the parametric variable 77 . In order to obtain this expression, notice that 


_ a = e (l/ 2 )ln(T- A) 

Hence, using equation (24) in this expression shows that for T e 



where 

mOj) s V cot 77 + A (26) 

It follows from the first equation (13) that the points in the physical plane (fig. 2) 
are related to the points in the T-plane by 


z(T) = 



dW Q 

- dT + Constant 

dT 


Substituting equations (19) and (23) into this formula shows that 


z(T) = — fl^T) + I 2 (T) 

77 L J 


(27) 


(28) 


where 


I X (T) 



T - 2 A + 2iVA VT - A 
T 



dT 


i 2 ( t ) 



T - 2 A + 2i Va Vt - A 


dT 


(29) 


(30) 
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and we have used the fact, indicated in figure 2, that the origin of the coordinate system 
in the physical plane is to be at the point D. On carrying out the indicated integration in 
equations (29) and (30) we find that 

Ij(T) = (T + 1) + 4iVA Vt - A - 4 A ln^iVT - A + Va^ 


4 Va Vl + A + 2Aln(Vl + A- Va^ + 4 A i?r 


(31) 


and that 

I 2 (T) = - 2iV^ V T „- a + 2 ln^iVT - A + V^) 


2 A + 2Va Vl + A - ln^VT + A - - 2?ri 


Hence, substituting these relations into equation (28) yields 


z(T) =i 

71 


_ 2 A 

T + —— + 


2iVZ^2 - l^VT - A + 2(1-2 A)ln^i V t - A + Va) 


(2 A + 1 ) + (2 A - l)ln(Vl + A - Va) 2 + 6Va Vl + A + 2i?r(2 A - 1) 


(32) 


By definition (see figs. 2 and 5), 


z(A) = a + ib 


Hence, equation (32) shows, after equating real and imaginary parts, 

- 2 

3(A + 1) + (2 A - l)ln ±A^ lA/A j + 6 Va VT 


a = — 


-| ^ 


+ A 


b = 2(2 A - 1) 


(33) 


(34) 
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There is an error in Ehrich's expression for a which corresponds to having the wrong 
coefficient in the Vl + A term in equation (34). By using these latter two equations, 
equation (32) can be written as 


z(T) = a + ib + 


(T - A)^l - + 2iVa^2 - IWt - 


A - b In 


lVt - A + VT 
Va / 


(35) 


We notice for future reference that equation (31) shows 


Ij(A) = A+ l + 2Aln + 4 Va Vl + A + 4?ri A 


(36) 


Formulation of First-Order Problem in Physical Plane 

The mapping T - z defined by equation (35) maps the upper half T-plane approxi- 
mately into the region of flow in the physical plane. The domain is mapped into the 
cross-hatched region of the physical plane shown in figure 6. The curve y q and the 
line segment EC are mapped into the dashed boundaries Sq and Fq, respectively, of 
this region. This region, of course, differs from the true interior of the jet, whose 
boundaries S and F are indicated by the solid lines (curved) in the figure. However, 
when |e| is sufficiently small, S and F do not differ very much from the zeroth-order 
boundaries of the jet S Q and Fq, respectively. 

Now the first group of the boundary conditions (11) are specified on the curves S 
and F in the physical plane, whose shapes are not known at this stage of the solution. 
However, the curves Sq and Fq, whose shapes are known, differ from S and F, re- 
spectively, by quantities which are of order e. Since insofar as the first-order solu- 
tion is concerned, the boundary conditions only have to be satisfied up to and including 



Figure 6. - Comparison of zeroth-order and true boundaries in physical plane. 
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terms of order €, we shall attempt to transfer the boundary conditions correct to terms 
of order e from S and F to Sq and Fq, respectively. To this end recall that the 
solutions have been divided into two parts (indicated by the superscripts + and -), one 
of which is holomorphic in D + and the other holomorphic in D~. We shall assume 
where necessary (as is done in thin-airfoil theory) that each of these portions of the 
solution can be analytically continued across S and F to Sq and Fq, respectively. 
Thus, the values of £* and W* at a point z® of S can be expressed in terms of their 

g 

values at the neighboring point Zq of Sq by performing a Taylor series expansion of 
these quantities about Zq. Of course, similar remarks apply to F. Hence, 



Substituting the asymptotic expansions (12) into these Taylor series and retaining the 
terms of 0(e) yields 


i 






w 




(37) 


(38) 


£ 






( 39 ) 
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W+ ( zF ) = w o( z o ) + £ [ w t( z o) + fo( z o) z f] + °( €2 ) < 40 > 

Hence, we have succeeded in obtaining expressions for the values of the dependent vari- 
ables W* and f* at the points of the unknown boundaries S and F in terms of their 
values on known boundaries Sq and Fq with an error of order e^. When these ex- 
pressions are substituted into the boundary conditions (11) and terms of O(e^) are ne- 
glected, expressions will be obtained which give the correct boundary conditions for the 
first-order solutions on the known boundaries Sq and Fq. It is clear from the form of 
these boundary conditions that it will first be necessary to obtain expressions for 
| g ± (z®) | ^ and | f ^z-F) | ^ correct to terms of O(e^) before this can be carried out. But 
it follows immediately from equation (37) that 
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Hence, 



In a completely analogous way, it follows that 



Substituting the expansions (38) and (40) to (42) and the last of expansions (12) into the 
first group of boundary conditions (11), equating the coefficients of e to the first power, 
and using the zeroth-order boundary conditions (17) to eliminate |?o( z ol| yield 



1 fM Zl F l=i 

(44) 

_?<K) ^o F r z K F J 2 

9m [wj(z®) + ”o( z o) z l] = 0 

(45) 

M*i«W4K]-o 

(46) 


(47) 
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Thus, equations (43) to (47) are the boundary conditions for the first-order solutions on 
the boundaries S and F "transferred" to the zeroth-order boundaries Sq and Fq. 
Hence, the first-order boundary value problem has been transformed from one in which 
the shapes of the boundaries are unknown to one in which they are known. Notice, how- 
ever, that these boundary conditions involve the variables W£\ z^, and z^ . But 

is completely determined in terms of by the second equation (13). In view of 
this the conditions (43), (45), and (48) may be thought of as two boundary conditions con- 
necting the variable with the variable £ j across Sq (or equivalently the variable 
with the variable W^) plus an equation which determines z® once ^ are known. 
Thus, subtracting equation (46) from (45) yields 

= J?#*W"(z®) (48) 


Then, in view of the second equation (13), equations (43) and (48) are the boundary con- 
ditions on S n which connect the solution ft in D + with the solution in D", and 
u S 1 + S 1 

equation (45) serves to determine once is known. (Actually, z^ will be deter- 
mined in a slightly different fashion. ) Similarly, the variable can be eliminated 
between the conditions (44) and (47) to yield a single boundary condition for (or W^) 
on Fq. To this end notice it follows from equation (13) that 


j_5o = jl_5o 

? 0 dz ? 0 dW 0 dz 


1_50 ? _ dln g 0 g 

5 0 dW Q 0 ' dW Q 


Hence, in particular 


But 




d0Q 

d( P o 


de o 
d( p 0 
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where 


6 n s arg C 


This shows that i(d In /dW Q ) p is real. (In fact, it shows that i(d In C 0 / dW o) F 

Z~—Z - 7. — V. 


is the curvature of the stream line Sq. ) Hence, 

Re 


1 „F 

- Re 

/ d ln ?o\ » / Z F\ F 

(-o')'* '~o F l J 


V dW 0 / z _ z F ^ ° j 
z - z o -1 


= Re 


(d In C 




dW, 


0 ; -o FL 


to «>r 


z-z r 


fd In C f 


= i 


dW, 


Re 


0 


■a-oV 


= i 



Using equation (47) to eliminate 


from this relation yields 
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Substituting this into equation (44) yields 








1 

2 


Or, dividing both sides by i(d In CqA* W o) F ’ 

z=z q 



is the boundary condition for £ ^ (or W^) on Fq. Thus, the boundary conditions on the 
curves Sq and Fq are given by equations (43), (48), and (49). The boundary conditions 
for the remaining (solid) boundaries are easily deduced by substituting the first asymp- 
totic expansion (12) into the second group of boundary conditions (11) and equating the 
coefficients of e to the first power. For convenience, we now collect in one place the 
complete set of first-order boundary conditions. 
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Transformation of First-Order Problem to Wg-Plane 

Now these boundary conditions completely specify a boundary value problem (or 
more precisely, two boundary value problems connected along the curve Sq) in the 
zeroth-order region of flow in the physical plane. However, under the change of vari- 
able z — W Q defined by equations (21) and (35), this boundary value problem can be 
transformed into one in the zeroth-order region of flow in the Wg-plane (which is indi- 
cated in fig. 3). The boundary conditions in the Wg-plane are 

n Rl< W 0> _ 1 

Jo< W 0> ? 0 (W oL 2 1 ? 0 (W Q 

^[wJ(W 0 ) - W“(W 0 )] = 0 
24 




W 0 e EC (52) 


Hm. L(W ) ^ - Wj(W )1 = ?(W) ^"1 - f 

L° ° <X 0 ? 0 (W 0 ) 1 °J [2 0 0 d! 0 J 


5 0 <W 0 ) 

A .^.0 

C 0 < W 0> 


W 0 GKD 


W Q e EK 


„ q(w 0 ) ^ 

9 m - - =0 W n GGD 

w 


(53) 


where the boundary conditions (17) have been combined with the boundary conditions (50) 
to obtain the boundary conditions (53). Naturally, the variables ^ and are defined 
on the unit strip below the real axis, and the variables and are defined on the 
upper half-plane. 

The two equations (51) can be combined to give a single jump condition for 
across DC. In order to see this, notice it follows from figure 3 and the definition of a 
derivative of a holomorphic function that, for Wq e DC. 

0 = -2- 9m[ WJ(W Q ) - W’(W 0 )] 

d( p o 


= 9m[w +(«p 0 ) - w-(«p 0 )] = ^[w^) - w- l( , o) ] 


= 9m 


dW, 


K(W 0 ) - Wi(W 0 jj 


Hence, using equations (13) shows that 
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0 = 9 m 


1 

~dWj 

C 0 < w 0> 

dz 


dWj 

dz 


= 9m 


gl(W 0 ) gj(W 0 )' 

Uo< w o> W 

Multiplying this result by i and adding it to the first equation (51) gives 

e}(w 0 ) ?j(w 0 ) x — 

1 U - U - 1 for W 0 S DC 


(54) 


W V W 0> 2 | C 0 (W Q ) | 2 

It is also necessary to transform the boundary condition (52). To this end notice that for 


W Q e EC 


dcp ( 


9m 


dW o ?T(w 0 ) - 

C 0 ( W o) J — 1 - - <(W 0 

d ^0 ?0 w o 


dcp. 


1 dW n 

2 u u d£ n 


Hence, 


9 m 


dcp. 


dW o ct(w 0 ) . 

C 0 (w 0 ) — - J — ^ - wtw 

° ° d? 0 C 0 ( W 0 ) ° 


= 9m 


i a dW n 

1 9 C 0 (w 0 ) 0 


2 dcp. 


d?r 


Therefore, using the definition of a derivative of a complex variable 


9m 


dW r 


dW o ?J(W 0 ) , 

UW Q ) — T - WiW 
° ° d? 0 ? 0 (W 0 ) 1 0 


= 9m — 


2 dW r 


M w o> 


dW„ 


dCr 


Differentiating by parts yields 


9m ^ ? 0 (W Q ) 


dW 0 ?i(W Q ) 


dW oL5o M w o>J 


^( w 0 ) dwf 


= 


1 + f 1 — ) 

_2 0 0 dW 0 \2 d? 0 /_ 


Upon using both equations (13) we find 
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-M?o (w o) — 


dW n 


dW 0 Cj(W 0 ) 


?t(w 0 ) 


dWj 


[d? 0 C 0 ( W 0)J M W 0> M W 0> dz 


lm ? n (W, „) -^-(i — Ufc(c 0 < w o> — 
0 0 dW 0 \2 d? 0 / \° 0 dW 0 


dW 0 ^(W 0 ) 

«0 V W 0>_ 


Therefore, 

j>= 0 for W Q e EC 

In view of this equation, let us define A + on the unit strip below the real axis in the 
Wg-plane by 


9m ? 0 (W Q )< 


_d_ ^0 
i dW 0 d ?0 


M _ i 

k 0 ( w 0 ) 2 _ 


A + (Wq) - ? 0 ( W 0 ) 


dWo 


dW o w w o) 


5o\W 2 /J 


-1 < t// Q < 0 


Clearly, A + is holomorphic on its domain of definition, and 


JUA + (W n ) = 0 for W 0 e EC 


(55) 


(56) 


Now it follows from boundary conditions (17) that ?q(Wq) is real for Wq e KD. Hence, 
it is readily seen from figure 3 and the definition of a derivative of a holomorphic func- 
tion that for Wq e KD 


dWQ 

d ? 0 


1 

(d? 0 /dW 0 ) 


1 

9u o 


where we have put 


u 0 = 

This shows that (dWQ /d^) is also real on KD. Hence, it follows from the first bound- 
ary condition (53) that for Wq e KD 
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dW. 


0 =■ 




d?o C 0 < W 0> d ?0 

Differentiating this result along HD 


^< W 0> 1 
5o< w o> 2 


= 9m 


f dW 0 

Cj(w 0 ) 

l«o 

jo< w o> 


o = 


;{'■ 

dW 0 


d ?0 

5o< w o> 

3 1 

dW 0 

5i< w 0> 

d( P o ( 

d ?0 

_<0< W 0> 

d 1 

dW 0 

?l< w o> 

dW 0 | 

d?o 

jo< w o> 


dcp. 


- 9m 


- 9m 


Therefore, in view of the fact that C Q (W 0 ) is real for W 0 e KD ’ we conclude from def_ 
inition (55) that 


9m A + (W 0 ) = 0 for W Q e KD 
A completely analogous argument shows that 

9 m A + (W q ) = 0 for W Q e KE 
Now define A on the upper half W^-plane by 


(57) 


(58) 


a “( w q) - S 0 < W 0> 


[ dW 0 

c ;< w 0 ) 

\^0 

? 0 < W 0 > 


- 2 


0<^ r 


(59) 


Then A - is holomorphic in the upper half-plane. An argument analogous to that used 
to deduce condition (57) suffices to show that 


9m A'(W 0 ) = 0 for W Q e GD 


(60) 
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It follows from the jump condition (54) and figure 3 that, for Wq e DC, 


?0< W 0> 


~dW 0 C|(W 0 )' 


dW ol5o ?o< w o>J 


<JW 0 q<w 0 ) 


dW 0 [ 5 o <0< W 0>J 


■ ?0< W 0> 


f dW 0 

q<w 0 ) 

i5o 

5o< w o> 


*?0<w 0 )i 3 


d( p 0 


dW n 


d? 0 |c 0 ( w o>! : 


Hence, in view of definitions (55) and (59) this shows that the jump condition of A* is 


A + (W 0 )- A“(W 0 )=i? 0 (W 0 )-i-<| 


dW„ 




< 0 (W 0 )| 2 C?(W 0 ) 


\ + 2 for W n 6 DC 


The right side of this equation can be put into a more convenient form as follows: Notice 
that, for W Q e DC, 


a < 

r 

dW Q 

1 1 

, - 3 < 

r 

1 dW 0 

1.1, 

d( P o 

d? 0 

_Ko (w o>l 2 ?o< w o>_ 

dcp Q 

<0< W 0> d <0 
V- 

? 0 (Wq) Wj 
J 


= -2i 


a 

l 

9 ^o 

?0< W 0> 


dW, 


Sji.4 


Hence, 


A + (W q ) - A“(W 0 ) = -i? 0 (W Q ) 


a 

l 

d( P o 

?o< w o> 


dW, 


0 9m — 


+ 2 


for W 0 e DC 


(61) 


and upon collecting equations (56) to (58) and equation (60) we have 
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(62) 


9m A + (Wq) = 0 


for< 


w o 

w n 


w„ 


e EC 
e KD 
eKE 


JWa~(Wq) = 0 for W 0 e GD 


(63) 


Hence, in view of equations (61) to (63), we have transformed the boundary value prob- 
lem in the Wp-plane for the holomorphic functions ^ (or W^) into a boundary value 
problem for the holomorphic functions A*. Clearly, once the functions A* are known, 
the original functions can be obtained by simple quadrature from definitions (55) and 
(59). However, it is necessary to transform the boundary value problem once more be- 
fore a solution can be obtained. 


Transformation of Boundary Value Problem for A ± into 
T-Plane and Solution of Problem 

Under the change of variable Wq — T defined by equation (21), the boundary value 
problem for A* in the Wq- plane posed by equations (61) to (63) is transformed into the 
following boundary value problem for the holomorphic functions A* in the T-plane 
(shown in fig. 5). 


a r dW n i 

A + (T) - A"(T) = -i? 0 (T) -2 ^ -±- 

0 ^oK ?0< T > 

u 1 

? 0 (T) 

+ 2 

9m A + (£ + i0) = 0 

i * -1 


9mh~m + i0) = 0 

T— ( 
1 

VI 


Clearly, the domains of definition of A + and A” 

are and 




(64) 


It is more convenient to work with the sectionally analytic function f2 defined on 
the upper half T-plane in terms of A"*" as follows 
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S2 + (T) s 


n(T) H 


dW. 


dT 


0 A + (T) 


O'(T) 


dW, 


dT 


° A"(T) 


>0 


(65) 


Since it is clear from equation (19) that dW^/dT is real when T is real, it follows 
from the conditions (64) that £2 must satisfy the following conditions: 


« (T) - «"(T) = T(T) for Te/ 


0 


9m £2(£ + iO) = 0 for -oo < £ < +°o 


( 66 ) 


where we have put 


T(T) = - 


dW n 


dT 


i? 0 ( T ) 


a 

1 

dcp Q 

5 0 (t) 


dW, 


^ 9m 


Cn( T ) 


for Tey, 


0 


(67) 


Now suppose © is a sectionally analytic function which satisfies the conditions (66). 
Then, if co is any function which is holomorphic in the interior of the upper half-plane 
and which is real on the real axis, the function 0 + co also satisfied the boundary condi- 
tions (66). First, suppose that co has no singularities on the real axis. Then, the 
Schwartz reflection principle shows that u > has an analytic continuation to the entire 
T-plane and therefore has a Taylor series expansion about the origin which has real co- 
efficients. Thus, co can be represented in the form 


m 


Z c n T ' 

n=0 


c n real 


If m = +°°, co has an essential singularity at T = +°°; if m is finite, <o has a pole of 
order m at T = +°°. It will be shown subsequently that the behavior of £ at T = °° 
dictates that m be finite. If co has singularities on the real axis, the requirement that 
co be real there shows that they cannot be branch points. Hence, these singularities 
must be poles or perhaps essential singularities. However, an investigation of the solu- 
tion shows that, if £ is to be bounded, the only singularity which can be allowed is a 
simple pole at the origin. In view of equation (19) the most general solution to the 
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boundary value problem (66) with certain restrictions of boundedness imposed can be 
written as 


m 


dW„ 


S2(T) = 0(T) + V c T n + a ^ Re T > 0 

/— / A ^jrri 


n=0 


dT 


( 68 ) 


where c n and a are real constants and 0 is a sectionally analytic function which 
vanishes at infinity, is bounded on the real axis, and satisfies the boundary value prob- 
lem (66). 

The function © can be constructed as follows: An investigation of the behavior of 
T(T) at T = °° and T = -1 shows that it vanishes at these points like some power of T. 
Hence, the Plemelj formulas (ref. 8) show that the Cauchy integral 

— f -^-dr (69) 

27ri t - T 
^0 

(where the integration is to be performed in a counterclockwise direction along y^) is a 
sectionally analytic function which is bounded on the real axis, vanishes at infinity, and 
which satisfies the jump condition (66). However, this function is not necessarily real 
for real values of T. But this can be compensated (as shown in ref. 8) by adding the 
function 



r(r) dr 
r - T 


to equation (69). (N otice that, if f is holomorphic in the upper half-plane, the function 
f defined by f (T) = f(T) is also holomorphic there, and f(£) + f (4) = f(ij) + f(£) is real). 
Thus, the function 0 with the required properties is defined by 


0(T) S 


1 j r(r)dr _ J_ / 

2iri / r - T 2 tu / 
"y r\ 


T(t) dT 
r - T 


ReT > 0 


(70) 


In order to complete the solution it is necessary to determine the real constants 
c n (n = 0, 1, 2, ... ) and a. Before this can be done, however, it is necessary to use 
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the preceding results to find as a function of T. In view of the mapping T — z de- 
fined by equation (35), which transforms the T-plane into the physical plane, this will 
complete the solution. 

Apropos of these remarks notice that definitions (55) and (65) combined with equa- 
tion (68) show that for T e @ + 


UT) 


dW, 


0 d 


dT dW„ 


dWp Aj(T) t 
d?n\?n(T) 2, 


oL s o 


m 


= 0(T) + £ c n Tn + 
n=0 


dW. 


a 


0 


dT 


(71) 


Since equations (19) and (20) show that dW Q /d? 0 = 0 at T = -1, integrating both sides 
of this equation yields 


dW n 




~^ T ) l]_ f T 1 

5q( t ) 2j ? 0 (T) 


0(T)dT 


m 


Ve n r T n d T + * r _l_ ^ 

/ y J_ x V T > J_ \ h™ dT 


dT (72) 


n=0 


Substituting in equations (19) and (20) and rearranging yield 


i^(T) i 


= — + 17T Va 


h™ 2 (T 


V T ) f T 1 
+ 1)Vt - A J_ x ? 0 (r 


— x — 0(T)dT 

? 0 ( T ) 


+ i 7 tVa 


? 0 ( T ) 


m 


(T + 1)Vt - A 

n=0 


• f — — T n dT 

J_ x M T > 


+ g ° (T> - F -L-!±±<n: (73) 

(T + 1) Vt - A J_-y M T > T 
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Since T(T) = O(T^) with some /3 < 0 for T — °°, it follows from the properties of the 
Cauchy -type integral (ref. 8) that 

0(T) = o(l) for T - * 


It is easy to see from equation (22) that £q(T) - 1 as T - <». Hence, as T - 

*T 


L 


L 


Cn(T) 


0(T) dT = o(T) 




1 rpn ^ T n+ 1 


and 


r 


1 T + 1 
Cn(T) T 


dT ~ T 


Using these results in equation (73) shows that 

m 


?J(T) 

^T) 





for T — 


Thus, in order that be bounded for T — it is necessary to set 

c n = 0 for n > 0 (74) 

Equation (73) can now be written as 


?i(t) 

MT) 


- + in Va 
2 


?n(T) 


(T + 1)VT - A 



0(T) dT + c Q 





1 T + 1 
? 0 ( T ) T 


dT 


? 0 (T) 

(T + 1)Vt - A 


X 



0(T) dT + c Q 



1 

? 0 (T) 


dT +■ 



1 

? 0 (T) 


T + 1 


dT 


T 


(75) 
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Since ?q(A) = -1 and since ©(A) is bounded, it follows from equation (75) that 


2(4 + 1) 


C 0 (T) 2 


£ 


?n(T) 


O(T) dT 


+ c, 


/•a r a 

/ _i_ dT + — / _1_ TjJL 

ii ? o< T) 'Jfj M T > T 


dT 


for T — A 


Hence, in order that £ ^ be bounded at T = A, it is necessary that Cq and a satisfy 
the two real equations 



O(T) dT + c Q 




1 

CoCD 


T + 1 
T 


dT = 0 


(76) 


Using equation (27) and the fact that the origin of the coordinate system in the physical 
plane is at D (fig. 2), it follows that 


z(T) 



1 

Co( T ) 


T + 1 


dT 


T 


(77) 


Hence, it follows from this and definition (29) that equation (76) can be written as 



0(T) dT + c 0 I 1 (A) + az(A) = 0 


(78) 


The constant oi appearing in this equation is directly related to a geometric parameter 
of the flow. In order to see this, notice first that the flow in the jet far downstream be- 
comes uniform. Hence, in view of equation (5) (see fig. 2), 

lim £ + (z) = |S + (z F )| = Vl + e 
z— c 
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Therefore, continuity requirements dictate that the dimensionless net volume flow 
through the jet Q be equal to 


nyr + e 

But Q must also be equal to the change in the stream function across the jet. Hence, 
the boundary conditions (11) show that 

o - = q = h yiTT 


Expanding both sides of this expression in powers of e by using the expansion (15), the 
last expansion (12), and equation (18) shows that 


1 - + • • • — (1 + ehi + . . .)(1 + — e+. , .) 

1 1 2 


Upon equating like powers of e, we find that 



(79) 


xp 

Another expression for if/^ can be obtained by first obtaining an expression for the 
first-order complex potential within the jet from equation (71). To this end, notice 
it follows from equations (13) that 


Hence, 


1 dW l dW l 1 dW 0 dW l 
? 0 Cq dz "dW Q 
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Substituting this result together with equation (74) into equation (71) yields 


dW 0 

”?j(T) _ i 

[din C 0 

C 0 (T) 2 


dW r 


Wt CT) + - Wq(T) } = 0(T) + c Q + a ” 

1 2 U I u dT 


Integrating this expression between -1 and T and recalling that dWp/d In Cq = 0 at 
T = -1 and that the normalization condition (7) implies Wj(T) = Wq(T) = 0 at T = -1, 
we find 


dW r 


dln? Q 


C"*"(T) ”1 f ,r ^~ 

— — - - - wt(T) +iw n (T) = / ©(T) dT + c n (T + 1) + aW n (T) (80) 

Uo< T > 2 J 2 h 


Using equation (72) together with equation (74) to eliminate £j(T)/£q(T) from this 
equation yields 


r— 


r /t 


r T 

|W+(T) -Iw 0 (T)J = 

1 1 0(T) dt — — j 

J x e 0 CT) c 0 (T) jf 

f 0(T) dT 

■1 


r r T 

— 


r 

z (T) - 1 

C 0 (T) 

+ c Q 

[l 

1 dT - 1 (T + 1) 

C 0 (T) C q (T) 

+ a 


(81) 


It follows from this equation, after using equations (76) and (77), that 
W+(A) = (I - «jw 0 (A) - ©(T) dT - c q (A + 1) 

Since 0(T) is real for real values of T (see fig. 5, p. 13) and since _^/«Wq(A) = -1, 
equation (82) shows that 

jWwt(A) = <* - - 

1 2 

An examination of figures 2, 5, and 6 (pp. 5, 13, and 17) together with equation (47) 
easily shows that 


(82) 


JWW+(A) = x^l 
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Hence, 


1 2 


Finally, combining this with equation (79) shows that 


a = -h. 


( 83 ) 


We now return to the problem of solving equation (78) for a and Cq, or, in view of 
equation (83), for and Cq. Substituting equations (23), (33), (34), (36), and (83) 
into equation (78) shows that 


/ 


T - 2A ± jliij gZI e(T) dT + c„ 

m « 


+ I + 24 * 4 yz VTT 

3(4 + X) + (24 - I)ln + 6 V 4 >T7" 


+ 4i7rAc 0 - 2ihj(2A - 1) = 0 


(84) 


Now 


/"A 

I T - 2A + 2i 

l 


2i Va Vt - A 


0(T) dT = 


/ A PL 

&U) d? - 2^A / 


(Va - iVT - A) 
T 


0(T) dT 


In view of the singularity in the denominator, the second integral must first be carried 
out over the path shown in figure 7, and then the limit 5—0 can be taken. After per- 
forming these operations, we find that 


2 Va J (V^ - iVt_-a) 0 ( T ) dT = 2 Va p. v. 


i:^ 




0(4) d4 - iff4A0(O) 


where the notation P. V. indicates that the Cauchy principal value of the singular integral 
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Figure 7. - Path of integration for 
in T-piane. 



(VA - i VT - A) 
T 


0(T)dT 


is to be taken. Substituting these relations into equation (84), recalling that 0(T) is 
real for real values of T, and taking the real and imaginary parts of the resulting ex- 
pression yield 



and 

27tACq - h^(2A - 1) + 2 ttA0(O) = 0 


Upon eliminating Cq between these two equations, we find 
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and 


c Q = (2A - hj - 0(0) (86) 

u 2 ttA 1 

Since equation (35) sets up a one-to-one correspondence between points of the physical 
plane and points of the T -plane, it is clear that equations (71) and (81) (with the use of 
definitions (67) and (70)) allow us to compute the first-order perturbation to the velocity 
and stream function at each point of the physical plane within the jet. Similar results 
are obtained by exactly the same procedure for the flow quantities in the main stream, 
but we shall not list the equations here. However, an inspection of these equations does 
show that the choice of the constants c n and a made above is sufficient to guarantee 
that is bounded. 

In view of the fact that once the shape of the jet is known it is quite easy to sketch 
in the streamline patterns, the most important quantities to be obtained from the analysis 
are the shapes of the curves S and F in the physical plane (see fig. 2, p. 5). IJow- 
ever, since the viscous spreading of the jet is controlled by the pressure (or equivalently 
the velocity) distribution along the slip line S, that quantity is also of some importance. 
Hence, explicit formulas will now be obtained for these quantities by using the formulas 
derived above. 


Derivation of Boundary Values 

In view of the one-to-one nature of the mapping involved, it is clear that, if 

W + (z) = i J/F, z must be a point on the stream line which passes through the point E 

in figure 2. In addition, since the velocity potential is increasing in the direction E— C 

along F, it is clear that, if /?eW + (z) > /?eW + (a + ib) (where a + ib are of course the 

coordinates of the point E), z must be a point of the free-stream line F. Hence, to 

o 

within an error of order e , the point z will be on the free-stream line F (i. e. , it 
will be on the first-order position of F) if 

^/«W + (z) = i//q + = -1 + 


and if 

/2<?W + (z) > /?eW 0 (a + ib) + e/?eW*(a + ib) 

where equation (18) has been used. It is also clear that /9eW 0 (zo ) > /?eW 0 (a + ib) for 
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TP 

any point zf e F n . In view of these considerations, it follows from equations (79), (82), 

uu / F \ FFF 

and (83) and the fact that 7 /«Wq(zq j = -1 that the point z = Zq + ez^ will be on the 

first-order position of F if z^ satisfies the equation 


W + (z F ) = W 0 (z F )(l - e^ F ) - 



0(T) dT - ec 0 (A + 1) 


(87) 


_L TP J. 

It is clear from this equation that, when Zq = a + ib, W (z ) = W (a + ib) and that 
$eW + (z F ) — as zf 1 — Hence, the point z F traverses the stream line F (to 

^ 2 TP 

within terms of order e ) as Zq traverses the zeroth-order free-stream line Fq. 

s s s 

Analogous considerations show that the point z = Zq + ezj, which is determined by 
the equation 


W(z S ) =W Q (z®)(l - (88) 

O 

traverses the first-order slip line S as the point Zq traverses the zeroth-order slip 
line Sq . 

By substituting equations (87) and (88) into the expansions (40) and (38), respectively, 

F S 

the first-order distances z ^ and z^ from the zeroth-order free-stream line to the 
free-stream line and from the zeroth-order slip line to the slip line, respectively, are 
found to be 


z i = -- 


Sn(T) 


W+(T) 


^ W 0 


(T) + 


i: 


O(T) dT + c n (A + 1) 


T e EC 


(89) 


and 


4 = L w i ( t) + ^ T e ^0 ( 90 > 

where the fact has been used that the curves Sq and Fq in the physical plane are the 
conformal images under the mapping T — z defined by equation (35) of the curve Sfo 
and the line EC, respectively, in the T -plane. The expansions (12) show that 

(91) 


z F = z(T) + ez F 


T e EC 
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z S = z(T) + ezj T e y Q 


o 

In addition, equation (41) shows that the magnitude of the velocity at each point z D 

2 

of the slip line is given to within terms of order e by 

q 9 9 9 [£t(T) 1 d£ n c 

| ?V) 1 2 = I C 0 cr) I : + 2 e 1 5 0 (T) I ■ Re -i- + -±- T e y Q (9 3) 

Cq(T) Sq(T) dz 


Upon substituting equations (89) and (90) into equations (91) and (92), respectively, and 
substituting equation (81) into the resulting expressions, we find, after using equa- 
tion (79), that 


z =z(T)(l + eh 




T ey n 


z F = z(A) + ^z(T) - z(A)J (1 + ehj) 
f /*T 


/ ^,h T) -„]«-^) flecrnc,,} 


dT> T e EC 


where equation (77) was also used in obtaining equation (95), 0 + (T) denotes the limiting 
value of 0(T) as T approaches y^ from within and the line integrals in equa- 

tion (94) are to be taken along y q 

It follows from the first equation (13) and equation (90) that 


?t(T) 

d?0 4 

5 0 (T) 

dW 0 

Cj(T) 

1 d ?o 

C 0 (T) 

C 0 (T) dW 0 

din C 0 

' <^0 A 

l 

o 

d In ^ \S 


[wJ(T) + ^[W 0 (T)J 


d In ? 0 

a Wq(T) 
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Upon substituting equation (80) into this expression, we find that 


Ct(T) 


d? 


os d _^of f 

z 1 dW o [J-1 


5 0 CT) Cq(t) dz uvv o 

Substituting this into equation (93) shows that 

|C + (z S )| 2 = \K 0 (. r T)\ 2 < 


0(T) dT + c q (T + 1) 


a 

2 


1+6 + 


2e/?e— f [e + (T) + c 0 ] 

<jw 0 J t L °J 


dT> T (96) 


where the line integral is along S^q. 

The distance S measured along the curve S is given by 


r 

s ^i 


dz k 


dT 


dT T e y Q 


(97) 


where the integral is taken along the curve SSq. Now, upon differentiating equation (94), 
we find after using the first equation (13) that 


dz 


S dW Q x 


dT dT 5 0 (T) 


1 + eh^ - e 


[e + (T) + c 0 ] 

dw o h L J 


dT> Tey, 


Hence, 


dz k 


dT 


?n(T) 


dW, 


dT 


f d in r T r . i 

h + 2eh t - 2e/?e ^ / 0 + (T) + c n 

1 dW 0 J.! 




dT 


1/2 


T ey n 


Substituting equation (96) into equation (79) gives, to within terms of order e , 


T e y n 


dz S 

= (l - 

dW Q 

1 

dT 

V 1 e ^i) 

dT 

|? + (z S )| 


Substituting this into equation (97) shows that 
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( 98 ) 


All the necessary results have now been obtained. However, it is convenient to 
rewrite some of these in more explicit form. This is done in appendix B. For con- 
venience, the most important equations of this section are now summarized. 


Summary of Equations 

Coordinates of downstream edge of nozzle: 


_ 1 
a TI 


3(A + 1) + (2 A - 1) In ^ V 1 ± A J 


+ A 


> 0 < A 


(34) 


b = 2(2 A - 1) 


Zeroth-order velocity: 


1 = T - 2A+ 2 i V a Vt - A () m T > o 


?n(T) 


(23) 


Transformation from T-plane to physical plane: 


;(T) + a + ib + i|(T - A) (l - + 2i Va^ 2 - i j Vt - A - b In ~ ^_ + 


(35) 


Change in stream function across jet: 


t// F = -1 - e/hj + h = 1- + ei// F 


(79) 


Asymptotic jet width: 


h = 1 + ehi 


(15) 
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where 


h i = 


- 27 T A 

(a+i)(4 a+i + ^-^a^aTi) 


2VIp.v. f* {Va±VaZl1 


©(4 )d| 


- ©( 4 )d 4 


A + 1 + 2 A In J + 4 VaVT 


+ A 


J 


0 ( 0 ) > ( 85 ) 


Intermediate variables: 


c = ( 2 A “ 1 ) . h. - 0(0) 
U 2 tt A 


( 86 ) 


1 T + 1 j 

1 J 

r tVt - a | 

i + lT + l\ 

77 T ' 

1< 0 (T) 1 

[Va (T + 1) 

v 2 T - A/ 


,(T) 


- (T + 1 ) Vt - A Re — — 


y T ) (t + i)Vt - aJ 


T = e -w? ( 99 ) 

sin rf 


Q + (v ) = - rfa) + P 


• v -- f 

Ziri / 

t/r\ 


r(y)M(y)dy 


27ri / y cot y - 77 cot 77 - i(y - 77) 

0 


27Ti / 

«/n 


r(y)M(y)dy 


y cot y - 77 cot 77 + i(y + 77) 


0 < 77 ^ 77 (B 4 ) 


where 


M(tj) 


sin^Tj 


(77 - cos 77 sin 77) + i 


(100) 
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T(y)M(y) 
y Cot y + ^ - iy 


dry 


-OO < £ < +°° 


zgfo) = z( — e" 11? 

\ sin rj 




•S 




_J_ e -1T ? 

sin rj 


J 


Vg(»)= l?V)|V«, 


J(tj) = (1 — 97 cot 77 + ir/) 



[ 1 — 7 ] COt 7] + A 


Position of slip line: 


z S = Zq(? 7)(1 + ehj) 



1 


0 + (y) + c qJ M(y)dy 


0 < ?] 


Velocity along slip line: 



? 0 <”> 



+ e 


+ 





0 < 7} < 7T 


(B5) 


(B6) 


(B8) 


(B9) 


(26) 


(B7) 


(BIO) 
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Distance along slip line: 


_S 

H, 


l 1 - € *r r v ~ 

o * J 0 Vg(y) 


V ,, , N 2 2 

°q (1 - y cot y) + y 


dy 0 <; 77 < 7 T 


(B13) 


y 


Position of free- stream line: 


z F = a + ib + [z(|) - (a + ib)] (1 + ehj) 


ru— l.' 

J ?o«l> ?o<*> 

A *- 


i) + c oK 


A £ ? < « (B14) 


Equations (19) and (20) have been combined with equation (Bl) to obtain equation (99). 


RESULTS AND DISCUSSION 

The numerical calculations were performed by using complex arithmetic. Hence, 
there is no need to separate the real and imaginary parts of the various formulas given 
in the preceding section. Equations (34) are used to calculate the orifice offset ratio 
B/A = b/a for various values of the parameter A. However, there is not a one-to-one 
correspondence between the values of B/A and the values of A. In view of this fact, it 
is more convenient to present the results in terms of the orifice orientation angle de- 
fined as tan - '*' B/A since there is a one-to-one correspondence between this latter 
quantity and the parameter A. A plot of the orifice angle against the parameter A is 
presented in figure 8. The orifice angle completely fixes the geometry of the problem. 
Hence, once the geometry of the orifice is set, the parameter A can be determined 
from figure 8. This parameter is the one which appears naturally in the formulas which 
are used to calculate the various physical quantities of interest. The only other param- 
eter appearing in the problem is e, which gives a measure of the difference between the 
total pressure in the jet and the total pressure in the main stream. This parameter is 
defined by equation (1) as 
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e 



Equations (23) and (113) are used to calculate T(ri) for various values of A, and 
these values of r(rj ) are used together with equation (114) to calculate 0 + (tj) and 0(£) 
for various values of A from equations (102) and (103), respectively. All the physical 
quantities presented in the plots are determined by these functions. Substituting the 
values of ©(£) into equation (85) allows the quantity h^ to be found as a function of the 
parameter A. Combining this with equation (15) gives h = H/Hq as a function of A 
and e. Now for two-dimensional jets, the jet contraction ratio is defined as the asymp- 
totic jet thickness divided by the length of the orifice. Hence, the jet contraction ratio 
is 


H _ h 

Va 2 + B 2 Va 2 + b 2 

Substituting equations (15) and (34) into this formula gives the jet contraction ratio as a 
function of A and e or, in view of figure 8, as a function of tan - '*' B/A and e. These 
results are presented in figure 9. It can be seen from figure 9 that, for positive values 
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of the orifice angle, small changes in e result in large changes in the contraction ratio, 
the effect becoming more marked as the orifice angle is increased. The opposite con- 
clusion holds for negative values of the orifice angle. Figure 9 also shows that, for a 
given orifice angle, increasing e always results in an increase in the jet contraction 
ratio. This increase is negligible, however, for orifice angles less than -100°. Fig- 
ure 9 shows that the jet contraction ratio is a maximum for an orifice angle of approxi- 
mately -80° and falls off markedly when the orifice angle is changed. 

The parametric equations (with parameter 17) for the slip line are obtained by sub- 
stituting equations (86), (35), (114), and (23) into equation (105) and using definitions 
(104) and the expression for 0 + (rj) discussed above. The parametric equations (with 
parameter £) of the free-stream line are obtained by substituting equations (23), (35), 
and (86) together with the values of 0(£) discussed above into equation (112). The re- 
sulting expressions (for the slip line and for the free-stream line) determine the bound- 
ary of the jet. The shapes of the jet boundaries for various values of the parameters e 
and B/A are shown in figures 10 to 19. A number of these figures consist of two parts. 
When this occurs, part (a) of the figure is for the case when the total pressure in the jet 
exceeds that in the main stream (e > 0) and part (b) is for the case when the total pres- 
sure in the main stream exceeds that in the jet (e < 0). Figure 10 corresponds to a jet 
injected normal to the main stream (B = 0). Figures 11 to 15 are for negative orifice 
angles (i. e. , jet injected downstream), and figures 16 to 19 are for positive orifice 
angles (i. e. , jet injected upstream). Figures 10 and 16 show that, when the orifice 
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6 , 




Dimensionless coordinate, X/Hq 
( b) Total pressure in main stream exceeds that in jet, e <0. 

Figure 10. - Jet contour for orifice offset ratio, B/A = 0. Jet injected normal to main stream. 


50 


Free-stream velocity, V c 



Figure 11. - Jet contour for orifice offset ratio, B/A = 2. Jet injected downstream. Total pressure in jet exceeds that in main stream, e > 0 


Free-stream velxity, 



Dimensionless coordinate, X/Hq 

Figure 12. - Jet contour for orifice offset ratio, B/A = -6. Jet injected downstream. 


Free-stream velxity, 



i I i i i i i 

LO 1.5 2.0 2.5 3.0 15 40 

Dimensionless cxrdinate, X/Hq 

Figure 13. - Jet contour for orifice offset ratio, B/A = -Z Jet injected downstream. 


4.5 


5.0 




Free-stream velocity, V e 



I i i i i i I i i 

0 .5 1.0 1.5 ^0 2. 5 3.0 3.5 4.0 

Dimensionless coordinate, X/Hq 

Figure 14. - Jet contour for orifice offset ratio, B/A - -1. Jet injected downstream. 



Dimensionless coordinate, X/H fl 
(bl Total pressure in main stream exceeds that in jet, « < 0. 

Figure 15. - Jet contour for orifice offset ratio, B/A * -0. 1862. Jet injected downstream. 





angle is greater than or equal to zero, a small change in the total pressure within the 
jet results in a large change in both the jet penetration and jet thickness. This effect 
becomes more pronounced as the orifice angle is increased. In fact, it becomes so pro- 
nounced that it is felt that the analysis may be invalid for orifice angles larger than 45°. 
For this reason, only the curves for e = 0 are given for these orifice angles. On the 
other hand, figures 11 to 15 show that, for sufficiently large negative values of the ori- 
fice angle, changes of the total pressure within the jet have almost no effect on the jet 
penetration or on the jet thickness. Although the scale of the figures is too small to 
show this, the numerical results do show that the jet always leaves the wall in a tan- 
gential direction. Figure 9 shows that turning the jet into the main stream tends to 
markedly decrease the contraction ratio, which is indicative of a decrease in the flow in 
the jet. The indications are, however, that a slight increase in the total pressure in the 
jet can easily compensate for this decreased flow. The jet penetration increases with 
increasing orifice angle. 

The pressure coefficient on the slip line is obtained as a function of the distance 
along the slip line in parametric form from equations (108) and (111) after using defini- 
tions (104), (106), and (107) and substituting in equations (114), (23), (79), and (86). 

The results of these calculations are shown in figures 20 to 27. Each figure is drawn 
for a different orifice angle. These curves contain all the information necessary for 
calculating the viscous boundary layer along the slip line. The figures show that the ve- 
locity within the jet at the upstream edge of the orifice decreases with both increasing 
orifice angle and decreasing e. It can also be seen from these curves that the pressure 
coefficient nearly reaches its asymptotic value in a distance of 10 jet diameters. 



Figure 20. - Pressure coefficient on slip line of jet for orifice offset ratio, B/A ■ 0. 
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Pressure coefficient, 










Figure 21 . - Pressure coefficient on slip line of jet for orifice offset ratio, B/A = 1 (first 
quadrant). Total pressure change within jet, e = 0. 

CONCLUDING REMARKS 

A procedure has been developed to obtain a solution to the problem of a two- 
dimensional jet injected from an orifice at an oblique angle to a moving stream for the 
case where the total pressure in the jet does not differ too much from the total pressure 
in the main stream. The analysis shows that, for orifices tilted into the main stream, 
small increases in the total pressure in the jet result in large increases in the jet pene- 
tration and in the jet thickness. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 24, 1969, 

129-01-07-07-22. 
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APPENDIX A 


FLUIDS OF UNEQUAL DENSITY 

A situation often encountered is the flow of a jet of high-density fluid into a low- 
density stream. When the jet and main stream are not composed of the same fluid, the 
foregoing analysis may be applied after a simple rescaling. This fact will now be proved. 
Previously, the problem was formulated in terms of the dimensionless velocity 


and its complex potential 


W = <p + iij/ = 


+ i'fr 
V ooH 0 


In order to consider fluids of unequal density let ^ be redefined as 


r 

Citi- flWriiY) ZSD+ 



Poo V oO 


S(z) = < 


if(z) = 


U - iV 
V 


z e d 


v_ 


(Al) 


where p ^ is the density of the fluid in the jet and is the density in the stream. To 
preserve the relation 


g _ dW 
dz 

it is also necessary to redefine the potential W as 
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W + (z) = 


z e D + 


W = < 


(# + i*k) 

P« V » H 0 


W~(z) = 

L 


$ + i'fr 

V H n 

oo 0 


z G D 


The dimensionless coordinates are, as before, 


x 



y 



a 



b 



h 



The perturbation parameter e is defined more generally as 


e 




= irV)! 2 - |e"(z s )| 2 


(A2) 


( 4 ) 


which shows that the jump condition on the slip line has the same form as for fluids of 
equal density. 

Likewise, the value of the dimensionless velocity on the free-stream line is given 
again by 
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(5) 








1 + € 


Very far from the orifice, the dimensionless velocity takes the same limits as in 
the equal density case, namely, 


S + (z) - 0 

rcz) - 1 

Similarly, the homogeneous conditions on the 
remain unaltered: 

9m ? + (z) = 0 

9m ? + (z) = 0 

9m £"(z) = 0 


z - K 

(9) 

z - G 

(10) 


velocity components at the solid boundaries 

z e KD 
z e EK 
z G GD 


Lastly, the boundary values of the stream function are again given by 


9mW + ( z) = 9mW~(z) = 0 z e S 
9mW + {z) = ij/ z e F 

•p 

The value of \p is to be determined in the course of the analysis exactly as before. 

Thus, it is seen that the boundary conditions on the velocities and the stream func- 
tion are identical to those obtained when the densities are equal (eq. (11)). Since the 
differential equation for the flow (Laplace's equation) remains unaltered as well, the 
analysis proceeds exactly as before, with the same solutions for the zeroth- and first- 
order problems. 

The interpretation of the results must be modified in one instance through inclusion 
of the density ratio. The pressure ratio (eq. (10)) is now 



p 0 ' p(x s> y s > 




KV )| 2 


(A3) 
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Equation (A3) shows that the pressure coefficient given in figures 20 to 27 must be re- 
interpreted as the ratio of the dynamic pressures when the fluid densities are unequal. 
Another quantity of interest is the mass flow in the jet 


“j = P ooVocHq 




i 


P f‘ 


(A 2 + B 2 )(l 


+ €) C.R. 


where C.R. is the contraction ratio computed to first order (fig. 9). 



APPENDIX B 


EXPLICIT FORMULAS FOR CALCULATING BOUNDARY VALUES 


Carrying out the differentiation by parts in equation (6'7) yields 


r(T) = - 


dWp 

dT 


C r- 


^ J_ + JL u(-L- JL . , 
?*o\So d V 5 0 Co d ^o Wo Wj 


T e y n 


It follows from the definition of a derivative of a holomorphic function that this can 
be written as 


F(T) = - 


dWg 

dT 




L dW oVo d V 5 0 <0 «0 \ dW 0 V 


- 2 




Tey 


0 

(Bl) 


Applying the Plemelj formulas (ref. 8) to equation (70), we see that 


0(T) - © + (T) = - T(T) + - P 


. V. f r(r)dr 1 / 

2ni / r - T 27ri / 


r ( T )dT _ _1_ / F(t) dr Tey 


T - T 


where the integration is to be performed in a counterclockwise direction along In 

view of equation (24), however, this can be written as 
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where, for brevity, we have written © + (77) in place of 0 + [(-77/sin 77) e _1? [], T(t 7) in 
place of r[-(T7/sin 77) e -iT [j, etc. Since 


d 




l y - cosy Sin y | 


sin 2 y 


equation (B 2 ) can be written as 



For T = £ + iO, equation ( 70 ) becomes 



and in view of equation (B 3 ) this can be written as 



©(£) = A 9**t 


s s 

Upon defining Zq(tj) and Zq(v) b Y 



y - COS y sin y . 


2 

Sin y 


dy 


y cot y + £ - iy 


< £ < + c 


(B5) 


e' 1 ”) 

Vsin n / 


ZH_ e - 1T J 


^sin rj 


) 


0 < r\ < 7r 


(B6) 


and using equation (B3), equation (94) becomes 


z S = z®fo)(l + chj) 


- e 


f V 1 

^ L'«m 


LC 0 (r) 5 0 (n)J 


[0 + (y) + c Q ]^ " C 0Sy Sin ^ + ijdy 0 < t? < 7r (B7) 
V sin*y / 


Upon defining VI{rj) by 

o 


V>) = V 00 |? + (z S )| o < T7 < 7 r 


(B8) 


and J(tj) by 


J(r]) = (1 - 7] COt 7] + ir/) 




(B9) 


where fj, is defined by equation (26), it follows from equations (19), (20), (24), (25), 
and (B6) that equation (96) can be written as 
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o<n<ir (bio) 


[v>)J 


K S 0 M 


1 + e 


v: 


+ 2ei t 


t e+M + c °]( 5 


rjl . cosy srny + 

. 2 
sm y 


Let p Q be the pressure far inside the orifice (the point K in fig. 2). Since the velocity 


is zero there, it is clear that p Q = P-. Hence, it follows from equation (2) that the 

s’ y s ) 


pressure at the point x y on the slip line p(x y ) is given by 

S S 


Po-P(yrs> _ |,+ (z |2 _ [ v sWT 

10 1 S I p 

± pV" v z 

^ oo OO 

Cl 


(Bll) 


Hence, let the pressure coefficient on the slip line C be defined by 

ir ^ 


C ps 


p 0 - p( x s > y s ) 

i * 1 - 


(B12) 


It is clear from equation (19) that 


dW 0 

= 1 

T + 1 

_ 1 ( | T | + 2 Re T + 1) 1//2 

dT 

77 

T 

7 T T 


Hence, it follows from equation (24) that 


dW 0 

_ sin n I 

( 2 

_JL_ 

dT 

ITT) i 

Vs in 2 


,1/2 


2r] cot 7] + 1 


for T Gy 


0 


Equations (24) and (101) show that for T e 
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dT = 


cot 77 2 i 


sin^77 


2 \V2 

dr] = — * — ( 1 - 2tj cot r] + —2 — J d?7 
shir l\ sin 2 r7, 


Hence, for T e Sf 


dW n 


dT 


0 


dT 


= — ( 2?7 cot 77 + 1 ] dTj = — [(77 cot 77 - l) 2 + 77 2 ld77 

^\sin 2 77 J n 7? 


Using this result in equation (98) together with definition (B8), 


(} - ^f) y . 


H„ 


r ■ 

y v; w 


|(1 - y cot y) 2 + y 2 


d y 


y 


Finally, it follows from equations (95) and (33) and figure 5 that 
z F = a + ib + [z(£) - (a + ib)] (1 + ehj) 

»«l 


- e f [— — l[o(4i) + c 0 ]d ?1 


A <! < 


(B13) 


(B14) 
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APPENDIX C 


A horizontal distance between edges 

of orifice 

a A/l 

B vertical distance between edges 

of orifice 

b B /l 

C _ pressure coefficient along slip line 
ps 

c fi n = 0, 1, 2, . . . constants 
D* flow regions in physical plane 

regions in T-plane 

F free streamline in physical plane 

H asymptotic jet width 

h H/l 

Ij function defined by eq. (29) 

J function defined by eq. (107) 

l characteristic length (set equal 

to H q ) 

M function defined by eq. (114) 

O order symbol 

o order symbol 

Pj total pressure in jet 

P ro total pressure in main stream 

P. V. Cauchy principal value 

p static pressure 

p Q static pressure at jet source 

p^ static pressure far upstream 

from jet 


dimensionless volume flow through 
jet 

S slip line in physical plane 

S distance along slip line 

slip line in T-plane 
T intermediate variable, i; + i r\ 

U X- component of velocity 

u U/V„ 

V Y-component of velocity 

Vg velocity along slip line inside of jet 
V M free-stream velocity 

" v/v„ 

W dimensionless complex potential, 

cp + 

X coordinate in physical plane 

x X/l 

Y coordinate in physical plane 

y y /i 

z dimensionless complex physical 

coordinate, x + iy 

g 

z dimensionless coordinate of points 

on slip line 
F 

z dimensionless coordinate of points 

on free-stream line 

ot constant 

T function defined by eq. (67) 

y dummy variable to replace rj 


SYMBOLS 
Q 
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A location of edge of orifice in 

T-plane 

« < p j - P J/f pvt 

£ dimensionless complex conjugate 

velocity, u - iv 

r\ coordinate in T-plane 

0 function defined by eq. (70) 

d o arg - ? 0 

A* defined by eqs. (55) and (59) 

Id function defined by eq. (2 6) 

£ coordinate in T-plane 

p density 

t dummy variable in T-plane 

$ velocity potential 

<P &/IV" 


if stream function 

* 

F 

^ value of dimensionless stream 

function on free-stream line 

£2 function defined by eq. (65) 

Subscripts: 

0 zeroth- order quantity 

1 first-order quantity 
Superscripts: 

F value of quantity on free- stream 
line 

S value of quantity on slip line 

+ value of quantity inside jet and 

orifice 

value of quantity in main stream 
complex conjugate (over bar) 
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